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ABSTRACT 

Combining archival data taken at radio and infrared wavelengths with state- 
of-the-art measurements at X-ray and gamma-ray energies, we assembled a broad- 
band spectral energy distribution (SED) of Cas A, a young supernova remnant. 
Except for strong thermal emission at infrared and X-ray wavelengths, the SED 
is dominated by non-thermal radiation. We attempted to model the non-thermal 
SED with a two-zone leptonic model which assumes that the radio emission is 
produced by electrons that are uniformly distributed throughout the remnant 
while the non-thermal X-ray emission is produced by electrons that are local- 
ized in regions near the forward shock. Synchrotron emission from the electrons 
can account for data from radio to X-ray wavelengths. Much of the GeV-TeV 
emission can also be explained by a combination of bremsstrahlung emission and 
inverse Compton scattering of infrared photons. However, the model cannot fit 
a distinct feature at GeV energies. This feature can be well accounted for by 
adding a 7r° emission component to the model, providing evidence for cosmic ray 
production in Cas A. We discuss the implications of these results. 

Subject headings: ISM: individual objects (Cassiopeia A), - ISM: supernova rem- 
nants, - radiation mechanisms: non-thermal 



1. Introduction 

Ever since their discovery, cosmic rays have been studied extensively. They consist of 
mostly protons and helium ions (but also heavy nuclei) and reach energies up to 10 21 eV. At 
energies below ~ 10 15 eV, cosmic rays are thought to be of Galactic origin. 



1 marayaar@purdue.edu 
2 cui@purdue.edu 



-2 - 



The most promising sources of Galactic cosmic rays are supernova remnants (SNRs). 
They are not only energetically favorable but also provide a mechanism to accelerate particles 
to very high energies. The process involves diffusion and scattering of particles across the 
shock front. The diffusive shock acceleration naturally leads to a power-law distribution of 
particles, which is observed in cosmic rays, and with an index that depends only on the 
compression factor of the shock (jBell 19781 ; iBlandford fc Eichler 19871 ). 



The link between SNRs and high-energy electrons was initially supported by radio ob- 
servations that showed synchrotron emission from these objects. The Chandra X-Ray Obser- 
vatory opened a new era of high-resolution X-ray imaging of SNRs. There is clear evidence 
for the existence of higher-energy electrons withi n these objects, wh i ch usually produce 



thin X-ray "filaments" associated with shocks (e.g.. iGotthelf et al. 2001r iHwang et al. 2002 



Long et al. 2003t iRho et al. 20021 ). The X-ray spectra of the filaments are consistent with 
synchrotron emission from particles as they move in an amplified magnetic field (e.g., 



Berezhko et al. 2002 ; IVink & Laming 2003 Berezhko fc Volk 2004 ; lArava et al. 2010h . 



Despite all the findings that link SNRs with high-energy electrons, direct evidence for 
the acceleration of ions by these objects is still lacking. In order to find a signature for 
hadronic acceleration, studying the photon spectrum of the sources at gamma-ray energies 
becomes fundamental. It is expected that at least a fraction of the emission at high en- 
ergies would be caused by the decay of neutral pions, which should be produced through 
collisions of accelerated hadrons with cold ambient protons and ions. Other mechanisms can 
also contribute to the gamma-ray spectrum besides hadronic interactions. The high-energy 
electrons that are responsible for the synchrotron emission at radio and X-ray wavelengths 
can upscatter ambient photons to TeV energies; and interactions between these electrons 
and other charged particles result in bremsstrahlung emission, which can also contribute 
to the emission in the GeV - TeV range. It then becomes important not only to obtain 
observations in the gamma-ray regime but also to model the broadband leptonic emission in 
detail, in order to separate and quantify the hadronic contribution. 

An example of a young SNR where high-quality broadband data are available is Cas A. 
The data should allow detailed studies of the properties of high-energy electrons in this rem- 
nant. For instance, radio images of Cas A indicate that the synchrotron emission is spread 
out within its shell, implying that the electrons responsible for it are likewise more or less uni- 
formly distributed. The X-ray emitting electrons are mainly confined near the forward sho ck, 
with their displacement being strongly limited by synchrotron loss (IGotthelf et al. 200l[ ). 



In this work, we present evidence for cosmic ray acceleration in Cas A from modeling 
its SED over a broad spectral range (from radio to gamma ray). 
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Data 



2.1. Archival data 



In order to construct a broadband SED, we used results from the literature, except for 
the GeV band where we analyzed public data from the LAT instrument on the Fermi satel- 
lite. The radio fluxes were me asured at the DKR-1000 radio telescope of the Pu shchino Ra- 
dio Astronomy Observatory (jVitkevich fc Kalachev 19651 ; lArtyukh et al. 1967) and at th e 
Staraya Pustyn' Radio Astronomy Observatory and taken from Vinyaikin (120061 ; 120071 ). 
at the frequencies 38, 151.5, 290, and 927 MHz for the epoch 2005.5. An additional 
Very Large Array flux at 74 MHz corresponding to the epoch 2005.2 was included 

( Helmboldt fc Kassim 2009 ). The VLA images of the remnant show that the radio 
is more diffuse, compared to the X-ray emission. 



emission 



The infrared fluxes were taken from t he Infrared Astronomical Satellite (IRAS) survey 
from 12 /im to 100 fim (IDwek et al. 19871 ). They are likely dominated by thermal emission 
from dust. Although we are only interested in non-thermal emission from Cas A, the infrared 
emission may be an important source of seed photons for leptonic gamma ray production. 

The X-ra y fluxes were obtain ed from deep exposures of Cas A with the Chandra X-ray 
Observatory (IHwang et al. 20041 ) . Data reduction was carried out as described in Araya 
et al. (2010). We excluded events from the central compact object in deriving the X- 
ray spectrum of the remnant. The X-ray spectrum of the remnant is dominat ed by thermal 
emission from highly ionized ejecta (IBorkowski et al. 19961 ; iFavata et al. 19971 ). as evidenced 
by the presence of numerous lines. Non-thermal X-ray emission manifests itself in thin 
filaments associated with the forward shock, as well as in the hard tail of the X-ray spectrum. 

The TeV gamma-ray measurements were taken from Acciari et al. (2010). They are 
based on observations of Cas A with VERITAS. 



2.2. Analysis of Fermi LAT data 



We analyzed the Fermi LAT data that were obtained b etwee n 2008 August 4 and 2010 
March 10. For details of the instrument, see Atwood et al. (120091 ) . 



1 The VLA array is an instrument of the National Radio Astronomy Observatory, a facility of the National 
Science Foundation operated under cooperative agreement by Associated Universities, Inc. 



'See the National Radio Astronomy Observatory image archive, http://images.nrao.edu/395 
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The data analysis was performed with the LAT Science Tools package §. The cuts 
and selection criteria applied to the data were those recommended by the Fermi Science 
Support Center. We selected diffuse events by excluding events with the lowest probability 
of being gamma rays with the gtselect tool and us ed the current vers ion of the instrument 



response functions given by PGJ/3-DIFFUSE ( jRando et al. 20091 ). The energy of the 
selected events ranges from 200 MeV to 300 GeV. We made appropriate time selections 
with the gtmktime tool. Specifically, we used a value of 105 for the zenith-angle cut to 
exclude times of high background that is associated with gamma rays produced in the Earth 
atmosphere. 

Due to the large LAT point-spread function at low energies, it becomes necessary to 
include events within a region around the source. This region of interest is taken here as 
a 10 radius circle around the cataloged position of Cas A. The accuracy of the analysis is 
improved when the contributions from other sources within a larger region ( "source region" ) 
are also included. We used a 20 ° radius source region, which is recommended for sources 
near the Galactic plane such as Cas A. The background included point sources at fixed 
positions locate d in the source reg ion, based on the recently released 11 month Fermi LAT 
Source Catalog jAbdo et al. 201(3 R as well as the galactic diffuse component (as specified in 



glLiem.v02.fit) and the isotropic extragalactic emission (as specified in isotropicJem-v02.txt) . 

Figured] shows a LAT image of the region near Cas A. A maximum likelihood method (as 
implemented in the gtlike tool) is usually used to determin e spectral paramete rs describing 



LAT sources of interest and their statistical significances (IMattox et al. 19961 ) . We found 
that Cas A was detected only at energies above 500 MeV. Between 500 MeV and 60 GeV, 
it was detected at a significance of 20.3 a. The best-fit position of the source, calculated 
with the gtfindsrc tool, was R.A. (J2000)= 23 h 23 m 21.9 s , decl. (J2000)= 58°49'47".4, with a 
statistical uncertainty circle of radius 0°.033 at the 90% confidence level. This position is in 
agre ement with that der ived by Abdo et al. (2010b). It also lies inside the VERITAS error 



box (lAcciari et al. 20101 ) 



The Fermi LAT spectrum of the source can be fitted with a power-law, dN/dE oc E p , 
with a photon index of p = 2.06 ± 0.07. The measured 0.5-60 GeV flux is (1.1 ± 0.1) x 10~ 8 
photons cm -2 s" 1 . Note that only statistical errors are shown here. The systematic errors 
on flux measurements have been estimated to be 5% at 500 MeV and 20% at 10 GeV 



( jRando et al. 20091 ). Our results are in overall agreement w ith those reported by the Fermi 



LAT team based on a smaller data set ( lAbdo et al. 2010bl ). 



3 See http://fermi.gsfc.nasa.gov/ssc 
4 See also http://fermi.gsfc.nasa.gov/ssc/data/access/lat/lyrjcatalog/ 
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3. Theoretical Modeling 

Figure [2] shows the broadband SED of Cas A. At infrared wavelengths, the SED is 
dominated by thermal emission from dust (which lies much above the extrapolated radio 
spectrum). Similarly, at X-ray wavelengths, the thermal emission from hot plasma stands 
out as the main component, as evidenced by the presence of numerous lines. In this work, 
we do not attempt to model the thermal emission, but instead focus on the radiation of 
non-thermal origin. 

At radio frequencies the spectrum of the source is of power-law shape, characteris- 
tic of optically thin synchrotron emission by high-energy electrons whose spectral energy 
distribution is also a power law. The existence of non-thermal electrons at even higher en- 
ergies has been most clearly reaffirmed by observations obtained with the Chandra X-ray 
Obse rvatory, which reveal the presence of thin compact regions of non-thermal X-ray emis- 



sion (IGotthelf et al. 20011 ). To provide a measure of the contribution of non-thermal emission 
to the overall X-ray spectrum, we also show, in Figure [21 the scaled-up X-ray spectrum of a 
representative non-thermal filament. Note that the X-ray spectrum is quite similar among 
non-thermal filaments near the forward shock (Araya et al. 2010). 

To account for the observed radio and X-ray emission, we constructed a model that 
involves two populations of electrons: diffuse electrons that are assumed to be uniformly dis- 
tributed over the entire SNR and localized electrons that are spatially concentrated near the 
forward shock. The synchrotron radiation from the former is expected to contribute mainly 
to radio and IR fluxes (but could also reach X-ray wavelengths if there are enough particles 
of sufficient energy) while that from the latter mainly to X-ray fluxes. Both populations 
may contribute to the observed gamma-ray fluxes via inverse Compton (IC) scattering and 
bremsstrahlung processes. 



3.1. Leptonic emission 

The two populations of electrons are modeled as zone 1 consisting of spherical "blobs" 
that approximate the non-thermal filaments seen, and zone 2 consisting of electrons uniformly 
distributed over a sphere of radius -Rsnr = 2-5 pc, assuming a distance <^SNR = 3.4 kpc 
faeed et al. 1995h . 

In reality, the radio emission shows spatial structures. Detailed spectral properties of 
the source such as a flattening of the spectrum above v ~ 30 GHz have been explained by 
Atoyan et al. (2000), who developed a two-zone model for the radio emission consisting 
of a diffuse component as well as localized "radio structures", such as an internal ring and 
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very compact knots. For simplicity, here we assumed that the synchrotron photons at radio 
frequencies are all associated with the diffuse electrons in zone 2. 

The description of zone 1 is based on our previous resul ts from a detailed study of 



non-thermal X-ray filaments near the forward shock of Cas A (lAraya et al. 20101 ) . In that 
work we used a synchrotron model that also takes into account the effects of diffusion and 
advection to explain the spectral evolution across the filaments. The analysis allowed us to 
determine the SED of the X-ray emitting electrons in the filaments as well as other properties 
such as magnetic field and level of turbulence. We found that the electron distributions are 
similar among the filaments examined. Here, we assumed that the electron populations in 
the blobs follow a power-law distribution of index 2.6, in a local magnetic field of 50 /xG. 
Note that the adopted value of the magnetic field is significantly lower than those estimated 
by as sociating the width o f filaments only with synchrotron cooling of radiating electrons 



(e.g. 



Vink fc Laming 20031 ). 



For zone 2, we modeled the electron distribution as a power law wi th a spectral inde x 



of 2.54, which reproduces the observed radio spectrum J u oc z/~ - 77 ( iBaars et al. 19771 ). 
The maximum energy of the electrons is poorly constrained by the data. To assess the 
maximum contribution from zone 2 to gamma-ray production in Cas A, we pushed the 
maximum electron energy to as high as possible without causing inconsistency with the X- 
ray measurement (see Figure [2]). A mean magnetic field of ~ 300 /zG was used to fit the 
data. This value is consistent with previous estimates derived from fits to the radio fluxes 
from the remnant (Atoyan et al. 2000). The overall normalization chosen gives an energy 
content in the particles that is also consistent with previous results (see Section 4). The high 
magnetic field could be caused by reverse shock (or second ary shocks) in the interior of the 



remnant (e.g., lLaming 2001t lUchiyama fc Aharonian 20081 ). 



Figure [2] also shows the results from leptonic modeling of the SED. The synchrotron 
emission from zones 1 and 2 is able to account for the non-thermal radio and X-ray fluxes. 
From the particle distributions, we calculated the number of gamma-ray photons emitted at 
the source per unit volume per second in the energy range e 7 to e 7 + de 7 as 

dn 7 (e 7 ) _ f°° <in 7 (T e ,e 7 ) ( ^\ / dJ . ^. 



dt J dt \ v e J \ dT 



e ■ 



where (4tt /v e )(dJ/dT) e is the number density per unit kinetic energy of electrons (T e ) and 
dn 7 (T e , e 7 ) / dt is the emissivity of one electron. The expected photon flux can then be 
determined for zones 1 and 2 (with volumes V\ and V2, respectively) as 



Vi 2 (2n 7 ' 2 (e T ) 



47 r ^SNR dt 



(2) 
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For the two zones, we included bremsstrahlung emission from electron-electron and 
electron-ion interactions, as well as IC scattering. The production rate of bremsstrahlung 
photons by an electron interacting with ambient protons, electrons, and helium ions with 
densities n p , n e , and n He , respectively, is given by 



cfa 7 (T e ,e 7 ) 
dt 



v e [(n p + An He ) a e!P (T e , e 7 ) + n e a ete (T e , e 7 )], (3) 



where a e e and the electron-electron and electron-ion cross sections, differential in 

photon energy and v e is the relative speed in bremsstrahlung collisions. The electron-ion 
cross section depends on the charge Z , which explains the factor n p + 4n# e in 

Equation 2. We considered only interactions of electrons with other electrons, protons, and 
fully ionized helium. This implies that n e = 1.2n p . We also assumed an abundance for 
helium of njje = 0.1 n p . 

The values of n p used in the calculations of expected bremsstrahlung fluxes can be 

(2) 

different for the two zones. For diffuse electrons in zone 2, the density n p is determined by 
the average density of protons in the remnant. If we consider uniform distribution of mass 
in a sphere of radius -RgNR,, then we get a mass for Cas A of 



M CasA~( 2 - 2 6M 



©J 



1 cm 



(2) 

Therefore, for an estimated remnant mass of ~ 10 M & (jWillingale et al. 2003[ ). we get n p ~ 
4.4 cm -3 . We fixed this value in the calculations. On the other hand, for the bremsstrahlung 
flux from electrons in zone 1 the density of ambient protons, n p l \ cannot exceed 0.8 cm -3 , 
in order for the model to be compatible with the Fermi LAT measurements (see Figure [2]). 
The fit shown in Figured was based on = 0.8 cm -3 . 

As for IC scattering, we considered various sources of seed photons. The most important 
one turns out to be infrared photons from the dust in the SNR. The IC scattering of these 
photons by electrons in zone 2 is calculated under the assumption that the infrared radiation 
is distributed uniformly across the SNR ( which is another gross simplification, as the IR 



intensity is seen to decrease going outward (IDwek et al. 19871 )). On the other hand, in order 



to fit the TeV fluxes, we found that about 15% of the total infrared photons reported by 
the IRAS observation should interact with the X-ray-emitting electrons located in zone 1 
near the boundary of Cas A. The IC scattering of other seed photons such as the cosmic 
microwave background (CMB) by electrons in zones 1 and 2, as well as the scattering of 
radio synchrotron photons from zone 2 by electrons in zone 1, were also included in the 
calculation. Figure |2] shows that these IC processes are only a minor contributor to the 
overall gamma-ray production in Cas A. 
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As also shown in Figure [21 even with the choice of extreme electron energies for zone 
2, the gamma-ray emission from this zone is still less than that from zone 1. It appears, 
therefore, that leptonic gamma-ray production in Cas A originates mainly near the forward 
shock. 



3.2. Need for a hadronic contribution 

It is clear from Figure [2]that our leptonic modeling fails to explain the observed gamma- 
ray emission at GeV energies. An additional component is needed to account for a distinct 
GeV feature of the SED. The presence of high-energy protons would naturally provide such 
a spectral component. Collisions between these protons and ambient (cold) protons lead to 
the production of pions. The decay of neutral pions produces gamma rays. 

Figure [3] shows the SED of Cas A with an additional component produced by 7r° decay. 
The g amma-ray spectrum from hadronic interacti ons was ca lculated as in Kamae et al. 



(120061 ). with a nuclear enhancement factor of 1.9 (IMori 20091 ) . The data seem to require 
that the spectral distribution of relativistic protons be a broken power law. In terms of the 
kinetic energy of the protons (T p ), the required distribution has the form dN/dT p oc T~ 2,1 
for T p < 11 GeV and dN/dT p oc T~ 2 ' 7 for T p > 11 GeV. 

The model shown in Figure [3J was obtained with the parameters in Table [H For con- 
sistency with the upper limit obtained by the Fermi LAT in the 200-500 MeV band, a 
local density of rip = 0.5 cm -3 was required for the bremsstrahlung emission from zone 1, 
which is lower than the value used for the model shown in Figure [2] (which has no hadronic 
contribution) . 



4. Discussion 



The Fermi satellite, launched on 2008 June 11, has opened a new window of exciting 
observations in the GeV regime. Hints of hadronic processes in SNRs have been seen. For 
instance, the Fermi LAT spect rum of the middl e-aged SNR W51C seems to indicate a 
dominant 7r°-decay contribution (jAbdo et al. 20091 ) . With respect to Cas A, the Fermi LAT 
team discarded the possibility that the GeV emission comes from t he compact object n ear the 
center of the remnant on the basis of spectral and timing analysis (jAbdo et al. 2010bl ). They 
carried out modeling of only the GeV-TeV gamma-ray spectrum in leptonic and hadronic 
scenarios separately. 



In this work, we modeled the radio and X-ray emission from Cas A as synchrotron radi- 
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ation from two populations of electrons. The radio emission is attributed to diffuse electrons 
distributed over the whole remnant, while the X-ray emission to electrons localized in about 
430 "blobs", each with a fixed diameter of 0.32 pc, which is a typical scale for the width of 
X-ray filaments near the forward shock. The combined volume of the "blobs" is about 1/10 
of the volume of the SNR. The total energy in the electrons is W e ~ 1.2 x 10 50 erK or abou t 
3% — 6% of the estimated kinetic energy in the explosion of Cas A (ILaming fc Hwang 20031 ). 



Bremsstrahlung radiation from the electrons can account for bulk of gamma-ray emission 
below about 1 GeV, while IC scattering of infrared photons by the electrons seems to be 
mainly responsible for gamma rays at TeV energies. Although there is considerable degen- 
eracy in modeling the radio to X-ray SED of the source, we found it difficult to account for 
the GeV spectrum of Cas A with leptons alone. 

Adding a 7r° component to the model can explain the distinct feature of the SED between 
about 1 and 40 GeV. The best match to the data is obtained for protons whose kinetic 
energies follow a broken power law with indices of ~ 2.1 for T p < 11 GeV and 2.7 for T p > 11 
GeV. Spectral steepening of the distribution toward higher energies is required to avoid 
exceeding the measured fluxes at TeV energies, where IC contribution is high. In our model, 
IC and 7r° components are comparable between 40 and 500 GeV, where the SED is relatively 
flat. The total energy in the protons is about W p ~ 0.6 x 10 50 erg. However, this value could 
get higher if a lower percentage of infrared emission is assumed to be interacting with the 
X-ray filaments (which would make IC contribution less important). 

Our model contains several simplifications. For instance, we neglected the spatial inho- 
mogeneities in the radio brightness distribution, which point at the existence of additional 
electron populations (Atoyan et al. 2000). However, as calculated by Atoyan et al. (2000), 
the expected bremsstrahlung flux from these electrons is considerably less than that of the 
diffuse component at GeV energies due to their steeper distributions. It can be argued sim- 
ilarly that the expected emission from IC scattering of ambient photons off these electrons 
is low, which justifies neglecting their contribution to the SED. We also assumed a uniform 
proton density inside Cas A, which is not realistic. We will investigate the effects of a more 
realistic density profile for young SNRs in a future work. 

Vink and Laming (2003) also modeled the emission from Cas A with two populations 
of electrons. They adopted a higher magnetic field (80-160 /xG) near the shock, which is 
obtained when comparing a typical width of an X-r ay filament with the synchrotron cooling 



length of a typical high-energy electron (see also iParizot et al. 20061 ). In our model we 
decided to use a lower value (50 /xG) derived by Araya et al. (2010) who took into account 
the effects of diffusion, in addition to those of synchrotron loss and advection. With their 
model, Vink and Laming (2003) predicted a gamma-ray flux from IC scattering of ambient 
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photons that lies below the HEGRA observation at TeV energies which, they argued, might 
be a hint to a pion decay origin of the radiation. In our model, most of the emission at 
TeV energies is attributed to IC scattering of infrared photons, which is a consequence of 
adopting a lower magnetic field at the shock front. It should be stressed that the data points 
at GeV energies, which were not available to previous studies, hold the key to revealing the 
presence of hadronic processes in Cas A. 

Although we believe that the lower value of the magnetic field that we adopted is more 
reliable, the measurement of the field is still quite uncertain. We also ran models with higher 
magnetic field. As an example, Figure 0] shows a reasonable fit to the data with B = 90 /zG 
in zone 1. The higher field causes a reduction in the IC contribution, which is compensated 
by a harder proton distribution (with a spectral index of 2.4 for kinetic energies T p > 8 
GeV). We also had to use a higher ambient proton density (n p = 1.5cm -3 ) and a lower 
volume for zone 1 (~ VsNR/^O)- For B > 90 /iG in zone 1, we found it difficult to explain 
the observed high-energy SED with our model. 

In summary, we were unable to account for the observed gamma-ray fluxes of Cas A 
at GeV energies with leptons alone in our model. The GeV "excess" can be well explained 
by additional gamma rays from the decay of neutral pions, which may be produced in the 
interactions between relativistic hadrons and dense matter in the remnant. This constitutes 
evidence for cosmic ray acceleration in Cas A. 

We thank T. Kamae for kindly providing us with his code for calculating the cross 
sections and gamma-ray fluxes in pp interactions, and S. Funk for providing his Fermi LAT 
fluxes of Cas A in electronic format for comparison. We also thank M. Lyutikov and D. 
Lomiashvili for useful discussions. This research has made use of data obtained from the 
Chandra Data Archive and the Chandra Source Catalog, and software provided by the 
Chandra X-ray Center (CXC) in the application package CIAO. This work has also made 
use of NASA's Astrophysical Data System. We gratefully acknowledge the financial support 
from NASA and Purdue University. 
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Table 1: Model Parameters 



Parameter 


Leptonic Component 


Hadronic Component 


Zone 1 


Zone 2 




Magnetic field (fiG) 


50 


300 




Ambient proton density (cm~ 3 ) 


0.5 


4.4 


4.4 


Minimum Lorentz factor 


1 


1 


1.5 


Maximum Lorentz factor 


7 x 10 7 


2 x 10 6 


5.5 x 10 5 


Volume 


~ T6 V SNR 


^SNR 


^SNR 


Spectral index* 


2.6 


2.54 


2.1, T p < 11 GeV 
2.7, T p > 11 GeV 


Total energy (erg) 


1.2 x 


10 50 


6.3 x 10 49 



-14- 




Galactic longitude 



Fig. 1. — Smoothed Fermi LAT count map (0.5-60 GeV) of a 16° x 16° region around Cas 
A. A binning of square pixels 0°.08 in size was used and smoothing was done (in ds9) with 
a Gaussian kernel of 2 pixel width. Contours are shown at the levels of 14, 28, 43, and 57 
counts. The positions of additional sources found in the Fermi Catalogue and included in 
the LAT data analysis are also indicated. 
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Fig. 2. — SED of Cas A and leptonic components. The dashed lines correspond to contribu- 
tions from zone 2 and the dash-dotted lines to those from zone 1. Shown in the GeV-TeV 
bands are bremsstrahlung fluxes (using a proton density for zone 1 of 0.8 cm -3 ), synchrotron 
self-Compton and IC scattering of CMB (IC/CMB), infrared (IC/IR), and radio (IC/radio) 
photons. For comparison with the model, we show, separately, the scaled-up X-ray spectrum 
of a non-thermal filament (in green). Note that the spectrum is "contaminated" by thermal 
emission at low energies (see Araya et al. 2010). 
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Fig. 3. — Modeling of the SED of Cas A with an additional hadronic component. The SED 
of protons is described by a broken power law in kinetic energy (see the text). For calculating 
the bremsstrahlung flux from zone 1, we adopted the value 0.5 cm~ 3 for the ambient proton 
density, as opposed to 0.8 cm -3 used in Figure |2j 
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